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Abstract 

A compact correlation-function expansion is developed for n'th order optical susceptibilities in 
the frequency domain using the Keldysh-Schwinger loop. By not keeping track of the relative time 
ordering of bra and ket interactions at the two branches of the loop, the resulting expressions contain 
only n + 1 basic terms, compared to the 2" terms required for a fully time-ordered density matrix 
description. Superoperator Green's function expressions for x^"^ derived using both expansions 
reflect different types of interferences between pathways .These are demonstrated for correlation- 
induced resonances in four wave mixing signals. 
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I. INTRODUCTION 

Time ordered expansions form the basis for the perturbative calculation of static and dy- 
namical properties of interacting many-body systems. The nonlinear response to a sequence 
of n short (impulsive) pulses is most naturally calculated in real (physical) time. The result- 
ing response functions contain 2" basic terms, stemming from the fact that each interaction 
can occur either with the ket or with the bra of the system density matrix. This fully time 
ordered expansion is routinely used for computing ultrafast (femtosecond) optical signals in 
molecules, semiconductors and other materials. The physical picture is recast in terms of 
the density matrix in Liouville space. Many-body theory of externally driven systems is in 
contrast commonly formulated using nonequilibrium Green's functions which act in Hilbert 
space [1-5]. Time-ordering is then maintained on an artificial Keldysh-Schwinger loop, [6,7] 
which corresponds to both forward and backward evolution in physical time and forms the 
basis for peturbative diagrammatic techniques. The loop provides a formal bookkeeping 
device for various interactions. We only keep track of the number of interactions with the 
ket and the bra but not of their relative time ordering. The nonlinear response function 
recast using these artificial (loop) time variables has then a considerably reduced number 
of terms, n + 1. Time-domain optical experiments performed using impulsive ultrashort 
pulses may be described on the loop, but the required transformation from loop-to real-time 
variables makes it hard to attribute physical meaning to the various terms [8] . 

In this paper we show that the loop time ordering is most suitable for computing nonlinear 
susceptibilities in the frequency domain, where real-time ordering is not maintained in any 
case. The frequency variables are directly conjugated to the various delay periods along 
the loop. In Sec. II we derive the correlation function loop expressions for the third order 
susceptibility. Since the loop expansion is much more compact, it may be advantageous 
to perform many-body calculations in the frequency domain on the loop and then switch 
to the time domain by a Fourier transform. This way one may exploit the full power of 
many body Green's function techniques. These expressions are then recast in Sec.III using a 
diagrammatic representation in terms of supcropcrators in Liouville space. The loop and the 
time-ordered expressions are compared in Sec. IV and shown to contain a different structure 
of resonances. A superficial look at the two types of expressions may suggest that they 
predict different types of resonances. This is however misleading since the various terms 
interfere. Gonsequently some apparent resonances may cancel and others may be induced by 
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dephasing processes. Simple diagrammatic rules are provided which allow to compute the 
partially time-ordered expressions. These subtle effects are illustrated in Sec.V by applying 
this formalism to study correlation-induced resonances in four wave mixing. The four point 
dipole correlation function is calculated for a multilevel system whose energy levels fluctuate 
by coupling to a Brownian oscillator bath. The model allows for an arbitrary degree of 
correlation between these fluctuations. The expressions may not be generally factorized 
into products of either real-time delays or loop-delays and the resulting complex pattern 
of resonances may not be attributed to specific time delays. When these fluctuations are 
negligible the loop expressions best reveal the resonances. In the limit of fast fluctuations 
(homogeneous dephasing) the real time expressions show these resonances. These subtle 
effects are demonstrated in Sec.V where we illustrate the different role of interference in the 
two types of expansion. We conclude by a discussion of these results in Sec. VI. 

II. NONLINEAR SUSCEPTIBILITIES ON THE KELDYSH LOOP 

We consider a system interacting with an external electric optical field E{t). The coupling 
Hamiltonian is Hint = —E{t)V, where V is the dipole operator. The nonlinear polarization 
P(") has n -Fl terms [9,10] 

n 

p(^)(t) = J2{i^^"~""\tW\i^^""\t))- (1) 

m=0 

Here | i/j^"^^ > is the perturbed wavefuntion to m'th order in the external field. We shall 
carry out the calculation for the third order response, n = 3. The generalization to n'th 
order is straightforward. The field consists of three modes and expanded as 

3 

E{t) = Y^ Ej{t) exp{-iujjt) + c.c. (2) 
i=i 

Eq. ([T]) now has four terms which correspond to m=3,2,l,0 and are represented by the 
Feynman diagrams (a), (b), (c), (d) shown in Fig. [1] respectively. The system interacts with 
the fields Ei, E2 and E3 at times ri, T2 and respectively, and the polarization is calculated 
at T4 by integrating over the time variables r,-.Each diagram represents a different ordering 
of Tj along the loop. 

Fourier transform of Eq. ([T]) to the frequency domain gives 

/oo 
dtexp{iuJst)P^''\t) = P^ + P, + P^ + P^ (3) 
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where 



/CO /•r4 l>T3 I-T2 

dn / dr^ / dT2 / dn E,{n)E2{T2)E^{Ts) (4) 



F{t4, T3, T2, Ti) exp(-ia;iTi - ia;2T2 - iuj^Ts + iuJsT^) 



du / dT2 / c^n / (ira E^{t^)E2{t2)E^{t^) (5) 
-00 J —00 J —00 J —00 



F{t3, T4, T2, Ti) exp(-ia;iri - ia;2r2 - iu^Ts + iuJsT4) 



/oo /■r4 /•T4 /•T2 

dT4 / dri / dT2 / dT3 £;i(ti)£;2(t2)£;3(t3) (6) 
-oo J —oo J —CO J —oo 



F{r3, T2, T4, Ti) exp(-ia;iTi - ia;2T2 - ia;3T3 + iujsr4) 

/oo l>T4 l>Tl pT2 

dn / dn / cir2 / rfrs Ei(ri)E2(T2)^3(T-3) (7) 
-oo J —oo J —oo J —oo 

F{tz, T2, Ti, T4) exp(-ia;iTi - ia;2r2 - ia;3r3 + iuJsTA). 
Here we have defined the correlation function 

F(t4,T3,T2,Ti) = {'^fTT[V{T4)V{n)V{T2)V{TM = C^f < V {t^W {ts)V {t2)V {t,) >, (8) 

where p is the equihbrium density matrix and V{t) are interaction picture operators with 
respect to the free system Hamiltonian H 

V{t) ^exp{iHt)Vex.p{-iHt). (9) 

By introducing Hcavyside step functions 6{t) we can set all time integration limits from — oo 
to oo, and combine the four terms as 

/oo /»oo /»oo /*oc 
/ / / dT4dTsdT2dTiexp{-iu;iTi - iu;2T2 - iu;sTs + iUgTi) (10) 
-OO J — OO J — oo J — oo 



-0{T2l)e{u2)0{U3)F{T3, n, T2, Ti )Ei (n) ^2 (r2)^3 (ts) 

+0{ni)d{n2)0{r2z)F{T^, r2, n)E,{n)E2ir2)EsiTs) 

-9{n^)9{n2)9{T23)F{Ts, T2, T,, n)E^{Ti)E2{T2)E^{n) 

The third order susceptibihty x*-^"* is defined by [9] 



P''-^\uJs) = j J J dUidU2dLJ3X^^\-Us]UJi,U2,UJ3)Ei{ui)E2{uJ2)E3{uJ3)6{uJs - UJi - UJ2 - UJ3), 

(11) 

where 

Ej{uj) = j dt Ej{t)exp{iujt). (12) 
By comparing Eqs. (fTOl) and (ITTl) we get 



2 /'oo roo roo 

{-Us]UJi,U2,UJ3) = ——^Y] / ^-^1 / ^-52 / C?S3 (13) 
l^^J p ^0 ^0 ^0 

[F{si + S2 + S3, Si + S2, Si, 0) exp(iCo'iSi + i{uJi + Co'2)s2 + iiuJi +UJ2 + CJ3)S3)- 

F{si + S2 - S3, Si + S2, Si, 0) exp(zco'iSi + i{uJi + 002)82 - - ujs + u;2)s3) + 
F{si - S2 - S3, Si - S2, Si, 0) exp(zu;iSi - i{uji - lUs)s2 - i{uJi - + co'2)s3)- 
F(0, S3, S3 + S2, S3 + S2 + Si) exp(iu;^si - ^(-a;^ + u;i)s2 - ii-ujs + uii + uj2)s^)] 

In Eq.(13) Sj are the time intervals between the various interactions along the loop. The 
frequency arguments of x^^^ are thus naturally connected with these variables. Time order- 
ing is thus maintained on the loop but not in real (physical) time, ^^denotes the sum over 
all 3! permutations of cui, ^2, i^3- 

We shall now compare this result with the fully time-ordered expressions for the response 
functions obtained by expanding the density matrix [8]. Generally S'''"^ has 2" terms. For 
n = 3 we get 
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oo 



p(3)(t) = / dti dt2 / dts (14) 
Jo Jo Jo 

S^^\t3, t2, ti)Ei{t -t3-t2- ti)E2{t -h- t2)E3it - tg). 



Unlike Eq. ( jlOl) . Ei, E2 and E3 now represent the first, the second, and the third pulse 
(chronologically ordered) 

s^'\h,t2,h) = C^fe{t,)eit2)eits) (15) 
< [[[v{h + t2 + h), v{t2 + h)], v{t3)], v{o)] > . 

Using Eq. ([HD, Eq. ([15]) gives 



^^'^^3, t2, h) = e{h)e{t2)e{h) [F(ti, + 12, h + t2 + o) (i6) 
+ F{o,h + t2,h + t2 + h,h) 

+ F{0,h,h+t2+h,ti+t2) 

+ F{U +t2 + ti + t2, ti, 0)] + c.c. 
The third order susceptibility is finally given by 

-j^ /'OO 1*00 poo 

X^^)(-uj,;uj^,uJ2,uJ3) = dts (17) 

S''^\t3,t2,ti) exp[iuJiti +i{uJi +0^2)^2 + i{uJi + uj2 + 003)13]. 

The loop expression Eq. f|T3|) generally has n + 1 terms (4 for n = 3) whereas the time 
ordered expression Eq. (fTTI) has a much larger number 2" (8 for n = 3). Note that the 
signal frequency does enter explicitly in the intergrations in Eq.(13) but not in Eq.(17). 
tj are intervals between successive interactions in real time and are most convenient for 
impulsive techniques, sj represent intervals along the loop and are particularly useful for 
frequency-domain susceptibilities. This will be demonstrated next. 

III. SUPEROPERATOR EXPRESSIONS FOR SUSCEPTIBILITIES 
By expressing Eq. ([8]) in terms of superoperators we can derive a more compact Green's 
function expressions for the susceptibilities using a simple diagrammatic representation. 
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Below we briefly survey the basic elements of the Liouville space superoperator formalism 
[8,11-13]. With each ordinary (Hilbert space) operator, Q, we associate two superoperators, 
denoted as Ql (left) and Qr (right) deflned through their left or right action on some Hilbert 
space operator X, 

QlX = QX, QrX = XQ. (18) 

We further define the linear combinations of these superoperators Q^={Ql + Qr) /2 and Q~ 
=Ql —Qr- Thus a +(— ) operation in Liouville space corresponds to an anticommutation 
(commutation) operation in Hilbert space, Q^X = {QX + XQ)/2 and Q^X =QX — XQ. 
The interaction picture for superoperators is defined by 

Uait) = exp{iLt)Ua exp{-iLt) a = L,R (19) 

where LA = [H, A] is the Liouville operator. 
When substituting Eq. (I9j) in Eq. ([H]) we get 

F{n,T3,T2,Ti) = i^fTr [eyip{-iHTu)V exp{-iHTi3)V exp{-iHr32)V exp{-iHT2i)V p] . 

(20) 

Note that all interactions in Eq. ([8]) are from the left i.e. they act on the ket of the density 
matrix. We can thus recast it using "left" superoperators as follows 

F{n, ra, T2, n) = C^fTr [\4(r4)\4(r3)FL(r2)\/L(ri)p] . (21) 

Combining Eqs. (fT9|) and fl20|) gives 

^(^4, r3,r2,ri) = 

{- fTr [exp {iLT4)VL exp ( -2^X43 ) Vl exp ( -iLr32 ) Vl exp ( -iLr2i ) Vl exp ( -iLri ) p] 

h 

When exp(— iLri) acts on the equilibrium density matrix p it does not affect it and 
gives p.Similarly exp(iLr4), when acts to the left will give 1 under the trace. These two 
propagators can thus be dropped and we finally get 

F{n,Ti,T2,n) = {^fTr [Viexp(-^Lr43)l^Lexp(-«Lr32)V^Lexp(-2Lr2i)rip] (22) 
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Using Eq. (1221) . the correlation functions in Eq. (fT3ll now become 

F{si + S2 + S3, Si + S2, Si, 0) =< VLg{s3)VLg{s2)VLg{si)VL > (23) 
F{si + S2- S3, Si + S2, Si, 0) =< VLg\s3)VLg{s2)VLg{si)VL > 
F{si -S2- S3, Si - S2, Si, 0) =< VLg\s3)VLg\s2)VLg{si)VL > 
F{0, S3, S3 + S2, S3 + S2 + Si) =< VLg\s3)VLg\s2)VLg\si)VL > • 

Here we have made use of the fact that all s variables are positive and represent "forward" 
propagation along the loop. 



g{s) = i-^Ms) exp(-zLs - rjs) (24) 

gHs) = (i:)^("5) exp(iLs — r]s) 
n 

We reiterate that ordering on the loop does not represent ordering in real time. Using 
superoperators we were able to recast F in terms of three propagators (in Hilbert space F, 
Eq. (12U]) has 4 propagators). By substituting Eq. (12^ in Eq. (IT^ we obtain 



< Vl^(u;i + CJ2 + U3)VLg{^^i + UJ2)VLg{uJi)VL > 

- < Vig^-UJs +OJI+ UJ2)VLg{uJl + U:2)VLg{uJl)VL > 

+ < VLg\-UOs + UJI+ UJ2)VLg\-UJs + UOl)VLg{uJi)VL > 

- < VLg\-uJs + cui + u:2)VLgK-u}s + Ui)VLg\-u:s)VL > 



Here 



gi^) = r— (26) 

UJ — L + tT] 

is the retarded Green's function, and 

g\^) = — ^ (27) 

uj — L — trj 

is the advanced Green's function, where t] is a positive infinitesimal. 
Eq.(25) may be represented by the unfolded loop diagrams shown in Fig. 2. These 
diagrams may be constructed using the following rules: 
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(i) Each Vl is represented by an arrow acting on the ket from the left. 

(ii) Each Vl is associated with one of the frequencies icui, ±co'2, icus ± ujg- Positive 
frequency +uj (negative frequency —u) is represented by an arrow pointing to the right 
(left). 

(iii) There are (n+ 1) choices for the position of Us along the loop see Eq. ([T]). Each 
gives one diagram. 

(iv) Each interval "before" ("after") Us gives a Green's function Q{uj), (^^(u;)). 

(v) The frequency argument of each Green's function is the sum of all "earlier" frequencies 
along the loop (frequency is cumulative). 

(vi) All ujj other than Ug can be interchanged, giving n! permutations oiuii-.-Un ■ Alto- 
gether x^^^ finally has {n + 1)\ terms. 

Finally, for comparison, using the fully time-ordered expansion Eq. f|T5|) we have 

x'^^\-u,-uJi,uj2,uj3) = (28) 

[{VLSiuJi +UJ2 + UJ3)VRg{uJi + UJ2)VRg{uJi)VL) 
+ {[VlG{uJi +UJ2 + UJ3)VrQ{uJi + UJ2)VlG{uJi)Vr) 
+ {VlQ{uJi +UJ2+ LU:i)VLQ{uJl + i02)VRQ{uJi)VR) 
+ {VlQ{uJi +UJ2+ U03)VlQ{uJi + UJ2)VlQMVl) 

- {VlG{uJi +U2 + U3)VlG{uJi + UJ2)VrQ{Ui)Vl) 

- {VlG{uJi +U2 + UJ3)VlG{uJi + UJ2)VlG{uJi)Vr) 

- {VlG{uJi +UJ2+ UJ3)VrQ{Ui + UJ2)VlG {uJi)Vl) 

- {VlG{uJi +U2+ UJ3,)VrQ{ui + UJ2)VrQ{ui)Vr)]. 

Eq.(28) may be represented by the double sided Feyman diagrams [8] shown in Fig.3. 
Note that this expression only contains retarded Green's functions representing forward 
time evolution, whereas Eq.(26) contains both retarded and advanced Green's functions. A 
more detailed comparison will be given in the next section. 

IV. RESONANCE STRUCTURE AND INTERFERENCE IN THE FULLY 
AND PARTIALLY TIME ORDERED EXPANSIONS 

Consider a multilevel system a, b, c... interacting with a bath whose Hamiltonian depends 
on the state of the system. The total eigenstates in the joint system -|- bath space are denoted 
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I aa >, I bp >, I C7 > etc. Note that the manifolds {a}, {P}, {7} diagonahze different bath 
Hamiltonians and they are not orthogonal to each other. For this model the total Green's 
function is given by 

gf^\ ^ I acy,bp »« aa,bp \ 

aa,bp ' 

where uJab = ^a — ^b is the transition frequency between states \a) and \b) and u^p = 6^ — 613. 
We next define the reduced bath Green's function by a partial trace over the system (denoted 
by a subscript s) 

g^^{uj) =« ab I g{uj) \ab», (30) 



Gabi^) is thus a superoperator in bath space 

I aP 

UJ - UJab- i^aP + If] 



^ \ap »« aP I , . 

^abMV — — (31) 

•i— I,]., I,] n-i-in 



a,l3 

Expanding Eq. (^5^ in eigenstates gives 

X^''\-UJs]UJ^,U2,UJ3) = --i- J2 P{(^)yadVdcVcbVba (32) 

a,b,c,d 

[ < Gdai^l + ^2 + UJ3)0cai^l + UJ2)QbaM > 
- < ^dai-^s +UJI+ U2)<3ca{'^l + UJ2)Qba{'^l) > 
+ < ^dai-^s +UJI+ UJ2)Gla{-UJs + ^^l)Gba{^^l) > 

Here P{a) is the equilibrium population of state \a). For comparison, by expanding the 
time-ordered expression Eq. fl25]) in eigenstates we get 

X(=^n-^s;^l,^2,^3) = -T^Yl E ^WKdK^cKbHa (33) 

p a,b,c,d 

[(^dc(^^l + <^2 + ^^3)Gdb{(^l + i^2)Gda{^^l)) 
+ {Qdc{^l +U}2+ U}3)Qdb{!^l + i^2)Qab{!^l)) 
+ {Qdc{^l +U}2+ U}3)Qac{^l + i^2)Qab{yJl)) 
+ {QbaiyJl + + ^^ZjQcaiyJl + ^2)'3da{.^l)) 

- {Qcbi.^1 + + ^^jQdbiyJl + ^2)Qda{.^l)) 

- {QdciyJl +UJ2 + UJs)Qdb{i^l + UJ2)Gda{uJl)) 

- {Qcbij^l + + ^3)Qca{!^l + ^2)Qda{MJl)) 

- {Qadiy^l +UJ2+ U}3)Qac{^l + ^2)Qab{i^l))] 
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It is interesting to note that Eq.(32) only suggests resonances with transitions involving 
the initial state a, uj^ai whereas Eq.(33) shows explicitly resonances between any pair of levels 
ojy^i. Both expressions are however formally exact and these apparent differences disappear 
by interference effects between various terms that can cancel some apparent resonances or 
induce new ones. When the dynamics of fluctuations is such that the loop time variables 
-51,52,33 are independent, the products of Green's functions in Eq. (123|) can be factorized 
and Eq. fl32|) then provides a natural representation for the observed resonances. Similarly, 
when the physical time delays between pulses tit2 and ^3 are independent, products of the 
corresponding Green's functions can be factorized, and Eq. ( 133|) should show the proper 
resonances. Most generally, neither factorization holds and averages of products of Green's 
functions must be carefully carried out. This will be illustrated in the next section using a 
model of multilevel system coupled to a Brownian oscillator bath. 

V. CORRELATION-INDUCED RESONANCES IN FOUR WAVE MIXING 

We consider a multilevel system coupled to a harmonic bath and described by the Hamil- 
tonian, 

H = Hs + Hb + Hsb (34) 
where the three terms represent respectively the system, the bath, and their interaction. 

Hs = J2e,\u){iy\, (35) 

V 

where Sy is the energy of eigenstate v. 

The system is linearly coupled to the bath through 

^SB = 5^4k)(i^|, (36) 

where Qy is a collective bath coordinate which modulates the energy of state v. 

The response function for this model of diagonal fluctuations can be calculated exactly 
using the second order cumulant expansion. Expanding the four point correlation function 
in the system eigenstates we get [14] 

F{t4, T3, T2, Ti) = (^)^ ^ VadVdcVcbVba eXp [-i{edT43 + ec1"32 + ebr2l) + fdcba{r4, ^3, ^2, Ti)] 

cba 

(37) 

where 
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fdcbain, T-3, Ta, Ti) = -^^^(^43) " ^cc(r32) - ^bfe(r2l) - gdcin2) 
+ Qdciji^) + 9dc{T32) - QdhijAl) + gdb{T42) 

+ Qdhim) - gdb{r32) - gcbirsi) + 5'cfe(i"32) + gcbim). 



(38) 



Here 

guu'it) = I dT2 I dTiC^„'{T2- n), (39) 
Jo Jo 

is the line broadening function, and 

a.'(r2-ri) = (g,(r2)g,.(ri)) (40) 

is the cross correlation function of frequency fluctuations of levels u and u'. We note the 
symmetry guu'{t) = gt'ui^^)- Upon the substitution of these results in Eqs. f|T3|) or f|T6l) and 
(fT7|) we can calculate x^^^ ■ Using the Brownian oscillator model for the correlations function 
we have [8] 

guu'it) = (^^^ - ^^)(exp(-A,,, |t| - 1 + A,,, |t|). (41) 

Here \^^, represents the coupling strength (the variances of frequency fluctuations are 
2\vyikT) and is the inverse timescale of bath fluctuations. For fast fluctuations >> 
2XkT we have 



9.At)=^..>\t\ (42) 
with r^^, = 2X^^,kT/A^^,. In the opposite limit of slow fluctuations A^ << 2XkT we get 

w(t) = A^^, kTf-tX^^,\t\ (43) 

Eq. (37) implies that Eq. fl23|l may not be generally factorized into three factors that depend 
on si, S2 and S3. Similarly Eq. (fT6|) may not be factorized into factors that depend on ti, t2, 
and ts; a three fold integration will be required to calculate x^^^ either representation. 

As an example for a dramatic interference effect related to these factorizations, let us 
consider the level system with a ground state a and two closely lying excited states b and d. b 
and d can represent, for example, two vibrational states belonging to the same electronically 
excited state or two Zeeman levels. The transition dipole only connects a with b and a with 
d. We look for two-photon resonances of the form (tui —UJ2— ^bd) in X^'^^ {^^s] —^2^^i)- 
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These are kind of Raman resonances but for excited state frequencies um- For simplicity 
we tune 003 to be off resonant and assume the fast fluctuation limit Eq. fj42|) . In this case 
the Green's functions assume the form 



The two diagrams responsible for such resonances in the time ordered expansion are 
shown in Fig. 4. They give the following contribution to x*^^"* 



,.1 1/1 1 

ab I I ^ ad 



Vab ?\ Vad P ^-.^ = -.^ ( + 



ixis - ^ha + i'^hh UI1-U2- Ulbd + i'^bd LJba + i^bb "1^2 " l^ad + i^dd 

(44) 



Here 



TM = {Ttb + rdd){l-v) (45) 
Abd = ArfrfAbb(l - r]) (46) 

7] is the correlation coefficient for fluctuations of and Uda- V = —1,0, 1 represent fully 
anticorrelated, uncorrelated and fully correlated fluctuations. The two terms in the brackets 
can be combined to give 

U1-UJ2- OObd + i{^bb + rrfd)(l -7] + 7]) 



(uJi - Uha + iTbb) - ^ad + i^dd) 

Substituting this in Eq. (H4l) results in 

1 1 

I T/ |2| T/ l2 
Vah Vnd 



(47) 



-UJs - OJba + i^bb OJl - UJba + i^bb -'^2 - '^ad + i'^ dd 
1 



1 + r7- 



(48) 



UJX-UJ2- ^bd + i'^bb + rrfrf)(l - r])_ 
The desired resonance (second term in the bracket) contains an 77 prefactor and its width 
scales as 1 — r]. When the fluctuations are uncorrelated (77 = 0) Eqs. (!23|) . (!38|) and (!25|) 
can be factorized and we expect no such resonances. Note that these resonances never 
show up in Eq. (|32l) but for 77 = cancel by interference in Eq. (|33ll . For flnite 77 and 
fast fluctuations Eq. ( l33l) can be factorized. Here, these resonances show up naturally in 
Eq. (l33l) . but in Eq. ( |32l) they come from the breakdown of the factorization. For 77 = 1 
the resonance width vanishes since there is no pure dephasing of the hd transition Eq. (l46ll . 
Such resonances have been observed both for coUisional broadening in atomic vapors [15] 
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and, for phonon broadening in mixed molecular crystals [16] and were denoted "dephasing 
induced" [15-17]. The present calculation shows that more precisely they are induced by 
the correlations of fluctuations rather than the fluctuations themselves, and are associated 
with specific factorizations of the multipoint correlation functions. 



When nonlinear response functions are calculated using the density matrix, the n'th order 
susceptibility has 2^n\ terms. These represent 2" Liouville space pathways which keep track 
of the complete time-ordering of the various interactions with the bra and the ket, combined 
with the n\ permutations of n frequencies representing all possible time-ordered interactions 
with the various fields. A wavefunction loop calculation keeps track of time ordering only 
partially (relative time ordering of ket and bra interactions is not maintained). This gives 
n + 1 terms, times the same n\ permutations for a total of (n+ 1)! terms. This considerable 
reduction in the number of terms is very convenient for the frequency-domain response, 
where the bookkeeping of time ordering is not necessary anyhow. 

When all field frequencies are tuned off resonance, we can neglect the imaginary part of 
the Green's function. We can then set Q=Q^ and use forward-only propagation. Eq. (!25|) 
then assumes a more symmetric form 



x('n-^s,^i,^2,c^3) = {VLGins)VLgin2)VLgin^)VL) sin, + n2 + ns + n^) (49) 



P4 denotes the summation over all 4! permutations of i72, ^3, ^4 with uji,uj2,uj3, —Us- 
Diagrams (a) (b) (c) and (d) correspond to the permutations —ujs = ^^4, f^s, and Qi 
respectively. Now we have a single basic term with [n + 1)\ permutations, as opposed 
Eq. fl2^ where we have [n + 1) terms each containing n\ permutations [18]. 

The loop expansion is most adequate for many-body perturbation theory and involves a 
combination of forward and backward time evolution periods in Hilbert space. The density 
matrix calculation, in contrast, only requires a forward propagation, but it must be done 
in Liouville space. The time-domain response functions may be obtained by a three-fold 
Fourier transform of x^^^ 



VI. 



DISCUSSION 




(50) 



X 



{~uJs]UJi,uj2,uJ:i)exp[-iuJi{ti + t2 + h) - iuj2{t2 + h) - iuj^h] 



14 



By substituting Eq. (l25l) in Eq. (150!) we can calculate the response function by transform- 
ing the compact frequency- domain expression obtained by a diagrammatic expansion on the 
Keldysh loop. 
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FIG. 1: The four loop diagrams for x*^'^^ representing Pa,Pb, Pc and P^ in Eq.(3). The loop expansion 
does not keep track of the relative time ordering of the bra and the ket. si, S2 and S3 are the time 
intervals ordered along the loop. 
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FIG. 2: Unfolded loop diagrams corresponding to diagrams (a), (b), (c) and (d) of Fig. 1. Eqs.(25) 
and (32) can be derived directly from these diagrams using the rules given in the text. All interac- 
tions are now from the left (ket), while the bra propagates freely. Solid and dashed lines represent 
forward and backward propagation, respectively. 
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FIG. 3: The eight double-sided Feynman diagrams representing the Liouville space pathways 
contributing to x^^^ Eq. (28). Complete time ordering of all interactions with the density matrix 
is maintained.il, t2 and ts are the physical time intervals between successive interactions, (a) and 
(d) of Fig.l are time ordered and each gives only one time ordered diagram (ai and di). (b) and 
(c) of Fig.l each split into 3 diagrams 61, 62 and 63 and ci, C2, C3. Altogether the four loop diagrams 
yield eight double-sided diagrams. 
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FIG. 4: The two double-sided diagrams which contribute to the correlation-induced resonance 
(Eq.44). 



19 



